Unit F3 
Integration 


1 The Riemann integral 


Introduction 


In this unit you will study the question: 


What do we mean by the area between the graph of a real function f 
and the x-axis? 


You will see how this can be answered by trapping the required area 
between increasingly accurate lower and upper estimates, each of which is 
the sum of the areas of suitably chosen rectangles. The area between the 
graph y = f(x) and the segment [a,b] of the x-axis is defined to be the 
supremum of the lower estimates and the infimum of the upper estimates, 
as long as these two values are equal. In this case, we call the common 
value the integral of f on [a,b], written 


fi or f toa 


You will see that, for many functions, we can evaluate integrals by using 
the Fundamental Theorem of Calculus, which allows us to think of 
integration as the inverse operation of differentiation. Although we will 
review techniques of integration, our main focus in this unit is on 
providing a rigorous foundation for the idea of integration, and on showing 
how this relates to concepts you have met in previous analysis units. 


Often it is not possible to evaluate an integral explicitly, and later in the 
unit you will meet methods for obtaining upper and lower bounds for the 
integral in such cases. You will also see how we can apply integration to 
derive some remarkable formulas for m and for estimating factorials, and to 
give a useful additional test for the convergence of certain series. 


1 The Riemann integral 


The purpose of this section is to give a rigorous definition of what we mean 
by the area between the graph 


y= f(x) (x€ [a,b]) 


and the closed interval [a,b] on the x-axis, and to explore its implications. 
We begin with an informal discussion to set the scene. 
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We have an intuitive notion of area. There are formulas for calculating the 
areas of simple geometric shapes such as rectangles and triangles, and we 
would certainly want our rigorous definition to agree with these. We would 
probably also agree that the region between the z-axis and the graph of a 
continuous function defined on a closed interval [a,b] has a definite area, 
even if we are uncertain how to measure it. On the other hand, for 
discontinuous functions it is not obvious that we can always say that the 
region between a graph and the z-axis has an area. For example, can we 
define such an area for the function 


zr’, 05741, 
2, 15142, 


f(x) = 


which has a discontinuity at the point x = 1? The graph of f is illustrated 
in Figure l(a). 
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(a) (b) 


Figure 1 The graphs of the discontinuous functions f and g 


As another example, can we define the area between the x-axis and the 
graph of the Dirichlet function on the closed interval [0,1], illustrated in 
Figure 1(b)? You met the Dirichlet function in Subsection 3.2 of 

Unit F1 Limits. Its rule is 


1, O0O<a< 1, z rational, 
g(x) = 


0, O<a2 <1, zx irrational, 


and in Unit F1 it was shown to be discontinuous at every point of its 
domain. 


It seems desirable that our definition of area should cover a wide range of 
functions. Later in this section we will prove that we can always assign a 
value to the area between the graph and the z-axis for a continuous 
function defined on a closed interval. You will also see that we can do the 
same for the function illustrated in Figure l(a), but not for the function 
illustrated in Figure 1(b). 


Our definition of area is based on finding lower and upper estimates for the 
‘area’ (if it exists) of the region between the graph y = f(x) and the 
x-axis, using the following approach. First we divide the interval fa, b] into 
a set of subintervals, called a partition of [a,b]. Then we consider two sets 
of rectangles, each rectangle having one of the subintervals as its base. In 
one set, we choose rectangles whose top edges lie on or below the graph, so 


the sum of their individual areas forms a lower estimate for the ‘area’ of 
the region; see Figure 2(a). In the other, we choose rectangles whose top 
edges lie on or above the graph, so the sum of their individual areas forms 
an upper estimate for the ‘area’; see Figure 2(b). 
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Figure 2 (a) A lower estimate and (b) an upper estimate for the area of the 
region between the graph y = f(x) and the z-axis 


In this way we can obtain many lower estimates and upper estimates by 
choosing different partitions of [a,b]. If there is a real number A with the 
properties 


the supremum of the lower estimates = A 
and 
the infimum of the upper estimates = A, 


then we define A to be the area between the graph and the x-axis. We call 
the number A the integral of f on [a,b], and denote it by 


[i or NECE 


We make all these ideas precise in the rest of this section, which is the 
longest and hardest section of the unit. On a first reading, you may wish 
to try to understand the main ideas without following every detail. The 
details may be easier to understand on a second reading. 


1.1 Definition of the integral 


In this subsection we work towards giving a rigorous definition of the area 
between the graph of a function f defined on a closed interval [a,b] and the 
x-axis; that is, the integral of f on [a,b]. 


Before we can give the definition, we need to introduce a number of key 
ideas. In the paragraphs below you will study: 


e some important terminology for functions 
e what is meant by a partition of a closed interval 


e the use of lower and upper Riemann sums to estimate areas. 
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Terminology for functions 


In Unit D1 Numbers you met the definitions of lower bound, greatest lower 
bound, upper bound and least upper bound of sets in R. Here you will 
meet analogous definitions for functions defined on an interval in R. These 
definitions of greatest lower bound and least upper bound generalise the 
notion of the minimum and mazimum of a function which you met in 
Unit D4 Continuity. The minimum and maximum of a function are 
illustrated in Figure 3 and we give a reminder of their definitions (in a 
slightly different form from those you saw in Unit D4) together with the 
definitions of lower and upper bounds for functions. 


YA YA 
M- 


(a) (b) 


Figure 3 (a) The minimum of a function (b) The maximum of a different 
function 


Definitions 


Let the function f be defined on the closed interval [a,b]. Then the 
following hold on fa, b]. 
e f is bounded below on [a,b] with m as a lower bound if 
f(z) =m, for all z € fa,b]. 
e mis the minimum of f on [a,b] if 
1. m is a lower bound for f on [a,b], and 
2. f(c) = m, for some c € [a,b]. 
Thus m = min{ f(x) : a <x < 6}, which we also write as min f or 
simply as min f. la; 
e f is bounded above on |a,b| with M as an upper bound if 
f(x) <M, forall z € [a,b]. 
e M is the maximum of f on [a,b] if 
1. M is an upper bound for f on [a,b], and 
2. f(d) = M, for some d € [a,b]. 
Thus M = max{f(x) : a <x < b}, which we also write as max f or 
simply as max f. [as 


e f is bounded on [a, }] if it is both bounded below and bounded 
above on |a, b|. 


1 The Riemann integral 


Note that any lower bound for a function f on [a,b] is less than or equal to 
any upper bound for f on [a,b]. 


A function f that is continuous on a closed interval [a,b] necessarily has 
both a minimum and a maximum, by the Extreme Value Theorem for 
continuous functions, as you saw in Subsection 3.3 of Unit D4. However, if 
f is not continuous on [a,b], then it may or may not have a minimum or a 
maximum on [a,b]; for example, the function in Figure 4(a) has neither a 
minimum nor a maximum on |a, b], whilst the function in Figure 4(b) has a 
maximum but no minimum on [a,b]. 


w= f(z) 


Figure 4 Two bounded functions on [a,b]: (a) a function with neither a 
minimum nor a maximum (b) a function with a maximum but no minimum 


Both functions in Figure 4 are certainly bounded on [a,b]; and there are 
numbers m and M, as shown, such that m < f(x) < M for all x € fa, b]. 
However, although the first function f takes values as close as we please to 
m and M, there is no point x in [a,b] where f(x) =m or M; and although 
the second function f takes values as close as we please to m, there is no 
point x in [a,b] where f(x) =m. This suggests the notions of greatest 
lower bound and least upper bound. 


Definitions 
Let the function f be defined on the closed interval [a,b]. Then 
e m is the infimum or greatest lower bound of f on {a, }] if 
1. m is a lower bound for f on [a,b], and 
2. if m >m, then f(c) = m’, for some c € [a,b]. 


Thus m = inf{ f(x) :a@< x < b}, which we also write as inf f or 


la, 


simply as inf f. 

e M is the supremum or least upper bound of f on [a,b] if 
1. M is an upper bound for f on |a, b|, and 
2. if M’ < M, then f(d) > M’, for some d € fa, b]. 


Thus M = sup{ f(x): a < x < b}, which we also write as sup f or 
[a,b] 


simply as sup f. 
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Remarks 


1. Note from the definition that the concepts of the infimum of a real 
function and the infimum of a set of real numbers are related in the 
following way: the infimum of a function is the infimum of the image set 
of the function. A similar remark applies to supremums. 


2. Any function f that is bounded on [a,b] necessarily possesses an 
infimum and a supremum on [a,b]: if m is any lower bound of f on 
a,b], then inf f > m and if M is any upper bound of f on [a,b], then 


sup f < M. 
3. If min f exists, then inf f = min f. Similarly, if max f exists, then 
[a,b] [a,b] [a,b] [a,b] 
sup f = max f. 
[a,b] [a,b] 


Worked Exercise F24 


Consider the function 


rozl =l] <z <1, 


a—a2, 15142. 


Sketch the graph of f and identify min f, max f, inf f and sup f (if they 
exist). 
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For each of the following functions f on [—1, 1], sketch the graph of f and 
identify min f, max f, inf f and sup f (if they exist). 


Partitions of a closed interval 


We now introduce the notion of a partition of a closed interval, which will 
play an important role in estimating the area between the graph of a 
function and the x-axis. 


1 The Riemann integral 
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Definitions 
A partition P of a closed interval [a,b] is a collection of a finite 
number of closed subintervals of [a, b], 
PN riper ile knalo a ki Pihos a ete 
where 
m= << a Koo Ky Ks KH, = 0. 
The points z;, 0£ i < n, are called the partition points of P. 


The ith subinterval is (1;—1,%;|, 1 < i < n, and its length is 
denoted by 62; = T; — 24-1. 


The mesh of P is the quantity || P|] = max KOZA 
<i<n 


A standard partition is a partition with subintervals of equal 
length. 


These definitions are illustrated in Figure 5. 


P 
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Figure 5 A partition P of an interval [a,b] 


Worked Exercise F25 


Let P be the partition of [0,1] given by 


P= {(0, 5], [5,3]. [5, a], [2 1 }- 
Find the mesh of P. 


1 The Riemann integral 


Exercise F34 


Write down the standard partition P of [—1,2] that contains 
four subintervals, and state the mesh of P. 


Lower and upper Riemann sums 


Next, we introduce the lower and upper Riemann sums for a bounded 
function f on an interval [a,b] with partition P; these correspond to 
underestimates and overestimates in our intuitive notion of the area 
between the graph y = f(x) and the z-axis. 


Definitions 


Let f be a bounded function on [a,b], and let P be the partition 
Lona Lans 2al 02 eee per Pomi Zol wie ng = 0 amd 
i, = lb. et 


mog = mi f) 390 E be seal 
and 

Mi supy fcr) oo E oi l 
for i= 1,2,...,n. 


Then the lower Riemann sum for f on [a,b] with partition P is 
n 
L(G) = eran 
i=1 
and the upper Riemann sum for f on [a,b] with partition P is 


OG = > Mada 
j= 


Note that the above definitions work equally well whether f takes positive 
or negative values in fa, b]. Regions between the graph y = f(x) and the 
x-axis where f takes negative values make a negative contribution to the 
lower and upper Riemann sums. In this subsection we will in general 
illustrate results for functions that are non-negative throughout [a,b], but 
we give further consideration to functions that are negative on all or part 
of the interval [a,b] in Subsection 1.3. 
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The terms m; and M; in the definitions denote the greatest lower bound 
and least upper bound of f on the ith subinterval of the partition; we need 
to use the infimum and supremum on the subintervals since the function f 
may not be continuous and so may not have a minimum or maximum on 
all (or any) subintervals. The lower Riemann sum is the sum of the areas 
of the rectangles with height m; and width ôx;, giving a lower estimate for 
the area between the graph and the z-axis, as illustrated in Figure 6(a). 
The upper Riemann sum is the sum of the areas of the rectangles with 
height M; and width 62;, giving an upper estimate for the area between 
the graph and the z-axis, as illustrated in Figure 6(b). 


YA YA 


y= f(x) 
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(a) (b) 


Figure 6 The rectangles whose areas are included in (a) the lower Riemann 
sum L(f, P) and (b) the upper Riemann sum U(f, P) 


Now, on any interval, the greatest lower bound of a function f is 
necessarily less than or equal to its least upper bound. It follows that, in 
each subinterval |x;—1, xi], we have m; < Mi. Summing from i = 1 to n, we 
obtain the following result. 


Theorem F43 


Let f be a bounded function on [a,b], and let P be a partition 
of [a,b]. Then 


L(f,P) < Uf, P). 


Worked Exercise F26 


Let 


20 O<a<il, 
f(z) = 
1, «=0,1, 


and let 
P= {(0, 4], [4,3], [> 4]> [451] } 
0, 


be a partition of [0,1]. Determine L(f, P) and U(f, P). 
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Solution 


®. It is helpful to make a sketch of the graph of f and the partition of 
the interval [0, 1]. 


0 T 3 
2 4 


í 
4 
Using the notation in the definitions, we set out the information 
needed to calculate the lower and upper Riemann sums using a layout 
that will make the calculation straightforward. Note that this 
function is increasing except at x = 0 and at x = 1, so we have to take 
special care at these points: the values of Mı and m4 are not what 
you might expect, and the values of mı and M4 are not taken by the 
function. © 


For the four subintervals in P, we have 


mı = 0, Mı = f(0) = ) 6x1 T 
m=fQ)=1, Mm=f@)=1, jm=h, 
m3 = IG = 1, M3 = f(#) = 5. 6x3 = Fe 
m=) Ma = 2, osa) 
Then 
4 
L(f,P) BE mol 
,= 1. 
= Org) a) a) xy) 
=O+}+3+4 
= 
8 
and 
4 
U(f,P) =>) Mid; 
j=l 
= (1x4) 


+ (1x a) + (2% 4) + (2% 3) 
3 
3 
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Exercise F35 


Let 

Seat 
and let 

P = {[0,3]; (5215 11} 


1 
5? 
be a partition of [0,1]. Determine L(f, P) and U(f, P). 


Worked Exercise F26 and Exercise F35 addressed the same function f on 
the same interval [0,1] but with different partitions. The two lower sums 
were 5/8 and 31/50, and the two upper sums were 11/8 and 3/2; each of 
the two lower sums is smaller than both of the two upper sums. 


We now look at another example to see whether this happens again in a 
different situation. We make use of the following formula for a sum of 
squares which you may have met in your previous studies (it can be proved 
by mathematical induction): 


1)(2 1 
12+22 +... +n? =>? = Mat Dene). for = l 2 ieee: 


Worked Exercise F27 


Let 


x, O0<¢< 1, 
2, 14142, 


and, for each n € N, let 
Pon = {[0, 5], [oah o 2- 2] 


be the standard partition of [0,2] into 2n equal subintervals. Determine 
L(f, Pan) and U(f, Pan). 
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Solution 
@. The graph of f is shown below. 


YA 
2 O, 
y = f(z) 
iA 
0 i SŠ 


The function f is increasing on [0,2]. Thus, on each subinterval in 
[0, 2], the infimum of f is the value of f at the left endpoint of the 
subinterval and the supremum of f is the value of f at the right 
endpoint of the subinterval. Since the ith subinterval in Pən is 


eer $ =Í 


= eee 
m=? ) m= s(4), =" = 
u n n n n 


Now f(z) = x? on [0,1] and f(z) = 2 on (1,2] so, since zp = I; we 
have 


2 
i=l . 
m= ( ) | fort = 12,5, NM 


i 
|, we have storm = Nh nea 2m 
n 


n 
and 
its =2, Word = FF 2, WFP Boa aay 2 


Hence 
n+1 


ID, 120) = Yom, = mu U M;0X; 


i=n+2 
= (i E) x + ((n—1) x2) x = 


eee 01 = 
£ (ei) 
n n n 


n(n + 1)(2n + 1) 2 
F S 
6n? n 
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Also, we have 


and 
M = 2 tor En n 2e 2m: 


Hence 


2n 
UGE Poy) — = M;62; 
za 


n 2n 
= So M;6x; + y MOT 
= 


i=n+1 


(ae oe T 
= (A ta 
n n 


= n(n + 1)(2n + 1) 
6n? 


We now look at the result of Worked Exercise F27 for different values of n. 
As n increases, the number, 2n, of subintervals in the partition increases 
and the length of each subinterval, 1/n, decreases. From the above 
formulas for L(f, Pən) and U(f, Pan), we find that, to three decimal places, 
the lower and upper Riemann sums are then as given in the following table. 


L(f, Pon) U(f, Pon) 
n=2: 4 equal subintervals 1.625 2.625 


n=: 8 equal subintervals 1.969 2.469 
n = 10: 20 equal subintervals 2.185 2.385 
n = 100: 200 equal subintervals 2.318 2.338 


As the subintervals increase in number and decrease in length, the lower 
sums increase and the upper sums decrease. But the lower sums are all less 
than or equal to all the upper sums! In fact this is always the case, as 
stated in the following result. 
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Theorem F44 


Let f be a bounded function on [a,b], and let P and P’ be partitions 
of [a,b]. Then 


ZU P) <U(f, P’). 


You saw in Theorem F43 that, for a given partition P, 
L(f, P) < U(f, P). 


In order to prove the general case in Theorem F44 with two different 
partitions P and P’ we require some new ideas. We develop these ideas 
and prove this result in Subsection 1.4. 


The integral 


We now return to our original problem: how to define an ‘integral’ that 
pins down our intuitive notion of ‘the area under a curve’. We have seen 
that lower Riemann sums provide underestimates for this ‘area’ and upper 
Riemann sums provide overestimates; see Figure 7. 
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YA 


y = f(z) 


> 
Zo Tı T2 T3 +++ Ln-1 Tn © 


Figure 7 Lower and upper Riemann sums 
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So we make the following definitions. 


Definitions 


Let f be a bounded function on a closed interval [a,b], let P be a 
partition of [a,b] and let L(f, P) and U(f, P), respectively, be the 
corresponding lower and upper Riemann sums. 


Then the lower integral of f on [a,b] is 


b 
E 
Ja 12 


and the upper integral of f on [a,b] is 


Fb 
f = inf U(f, P). 


We say that f is integrable on [a,b] if 


[ss fo 


and their common value is then called the integral of f on [a,b]. 


b b 
The integral is written as | f or i f(a) dx, and a and b are called 


the limits of integration. 


Remarks 


1. The common value of the lower and upper integrals of f (when it exists) 
is sometimes known as the Riemann integral of f, rather than simply as 
the integral of f. 


2. In the definitions, sup L(f, P) is the supremum of the lower Riemann 
P 


sums over all possible partitions P of the interval [a,b]. Similarly, 
inf U(f, P) is the infimum of the upper Riemann sums over all possible 


partitions P. 


3. For any bounded function f on an interval [a,b], it follows from the fact 
that all lower Riemann sums are less than or equal to all upper 
Riemann sums (Theorem F44), and from the above definitions, that the 


b b 
lower integral f f and the upper integral / f both exist, and that 
a a 


we always have 


[ts ft 


(We omit the proofs of these facts.) Note, however, that the lower and 
upper integrals take different values unless f is integrable. 
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We now return to the function f that we considered earlier in Worked 
Exercise F27. 


Worked Exercise F28 


Let 


2 
Prove that f is integrable on [0,2], and evaluate f f. 
0 


Solution 
®. It is helpful to sketch the graph. 


YA 

2 oeo 
y = f(x) 

la 


We have already seen in Worked Exercise F27 that if we take Pən to 
be the partition 


ET MJ) 


of (0, 2], then 
n(n + 1)(2n + 1) 2 
and 
n(n + 1)(2n + 1) 
U(f, R) = 6n3 2 
Then, as n + co, we have 
1+1 2+1 2 1 
JV oa oes 
6 mo g 3 


so that, in particular, 


[ret 
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2 
®. This holds because l f is defined to be the supremum of the 
0 


lower Riemann sums of f over all possible partitions of [0, 2], so it 
must be greater than or equal to the particular lower Riemann 
sum L(f, Pan) for any value of n. 8 


Similarly, as n — oo, we have 


Ci au MORE. Pa eel 


so that 
T 
[rss 
0 3 
We have now shown that 
2 aD 
te i f< i ged 
3— Jo 0 3 
It follows that 
2 
a 
Jo 0 


2 
so f is integrable on [0,2] and a J= L 
0 


However, not all bounded functions defined on closed intervals are 
integrable! 


Worked Exercise F29 


Let f be the Dirichlet function on [0,1] defined by 


1, O<a<1, x rational 
f(x) = are 
0, O<a <1, z irrational. 


1 71 
Determine the values of f f and J f, and hence show that f is not 
Jo 0 


integrable on [0, 1]. 
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Solution 


®. The graph of f is shown below. 


YA 
] Pcceccccccccccccccccccccccce e 
—Oeseseesesccscesssecescessse > ve 


Let P= they a wie tole ewe). herero — 0; a7-— 1 be 
any partition of [0,1]. Then, on each subinterval (x;—1, x;| in P, we 
have 

file =O ail Mb il, nor = 1,2 coon ft 


®. This is because every subinterval contains both rational and 
irrational points, by the density property of the real numbers (see 
Subsection 1.4 of Unit D1). @ 


So 


LF P) = MjOX; = 5 (0 x Om) = |) 
zli 


i= 
and 
COAR) = z Mi;ôzxi = ) (1 x O) = 2 = il 
j=l j= g= 


since the sum of the lengths of all the subintervals is equal to the 
length of the interval [0, 1]. 


It follows that 
1 
f f=suw1(f,P)=0 
Jo P 
and 


71 
i f =inf U(f, P) =p 


1 71 
Then, since | TFs f f, we conclude that f is not integrable 
JO 0 


on [0,1]. 
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Exercise F36 
For the function 
f(z)=12, xe [01], 
and the standard partition of [0, 1] 
r= {otl "afb nen, 
n| |n n n 


determine the values of L(f, Pa) and U(f, Pa). Hence show whether f is 


1 
integrable on [0,1] and, if it is, determine the value of f f. 
0 


1.2 Criteria for integrability 


It would be tedious to have to go back to the definition of integrability 
whenever we wish to show that a given function is integrable on a closed 
interval. You will now meet a number of criteria that we can use to avoid 
this. 


In order to prove that a bounded function f is integrable on [a,b] directly 
from the definition of the integral, we need to look at sup L(f, P) and 
P 


inf U(f, P) over all partitions P of [a,b]. However, as you saw in Worked 


Exercise F28 and Exercise F36, in many situations it is sufficient to 
consider just one sequence of partitions (P,,) in order to establish 
integrability. The conditions under which this simplification holds are set 
out in the following result. (Recall that ||P,,||, the mesh of Pa, is the length 
of the longest subinterval of Pn.) 


Theorem F45 


Let f be a bounded function on [a,b]. If there is a sequence of 
partitions (P,,) of [a,b] such that ||P,|| > 0 and 


lim L(f,P,) = lim U(f,P,) =A, where AER, 
noo n—> oo 


b 
then f is integrable on [a,b] and / j = A 


Proof Let e > 0. It follows from the equations in the statement of the 
theorem that there exists an integer n such that 


L(f,Pa)> A— łe and U(f, Pa) < A+ 4e. (1) 


Now, by the definitions of upper and lower integrals, 


b Fb 
f EIHP, and Ee (2) 


1 The Riemann integral 


Combining inequalities (1) and (2), we obtain 
b Fb 
A- łe< f f< f fSA+je. 


®. Here we have used the inequality [ JE [ f mentioned in the 


remarks after the definition of the integral in the previous subsection. As 
stated there, this follows from Theorem F44. 8 


Since £ is any positive number, we deduce that the upper and lower 
integrals of f on [a,b] are equal to A, so f is integrable and 


[tra a 


In fact the following result in the opposite direction to Theorem F45 also 
holds. However, its proof is somewhat more complicated, and we defer this 
to Subsection 1.4. 


Theorem F46 


If f is an integrable function on [a,b] and (Pp) is a sequence of 
partitions of [a,b] such that ||P,|| > 0, then 


dim Lf, Pa) = lim U P a Jo 


This result is particularly useful for proving that a function is not 
integrable. For if f is defined on [a,b] and we can find a sequence of 
partitions (P,,) of [a,b] whose mesh tends to zero but for which 


lim L(f, Pa) £ lim U(f, Pa), 


then it follows from Theorem F46 that f is not integrable on [a,b]. 


Exercise F37 


For each of the following functions f, determine whether f is integrable 
1 


on [0,1] and, if it is, find f f. 
0 


=). 0<z<1, 
a t)= 
(a) Fa) R nm 
b) f(x) xz, O<a<1, zx rational, 
4 64 = 
0, O<a<1, r irrational. 


141 


Unit F3 


142 


Integration 


The next result is of particular interest in that its statement says nothing 
about the value of the integral itself: it mentions only the difference 


between the lower and the upper Riemann sums. The result follows from 
Theorems F45 and F46. 


Corollary F47 Riemann’s Criterion 
Let f be bounded on [a,b]. Then 

f is integrable on [a, b] 
if and only if 


there is a sequence (P,,) of partitions of [a,b] with ||P,|| + 0 
such that U(f, Pa) — L(f, Pn) > 0. 


Proof Theorem F46 implies that if f is integrable on [a,b], and (P,) is a 
sequence of partitions of [a,b] with ||P,,|| > 0, then 


On the other hand, if there is a sequence (P,,) of partitions of [a,b] with 
||P, || + 0 such that U(f, P,) — L(f, P,) + 0, then, because 


L(f, Pa) zli < fr < U(f, Pa); 


these upper and lower integrals must be equal, so f is integrable on [a,b], 
by Theorem F45. E 


Two classes of integrable functions 


With Riemann’s Criterion at our disposal, we can now determine some 
large classes of functions that are always integrable: the monotonic 
functions and the continuous functions. 


Theorem F48 


A function f which is bounded and monotonic on fa, b] is integrable 
on [a,b]. 


Proof We prove this theorem in the case when f is increasing on fa, b]. 
(The proof is similar if f is decreasing.) 


Consider the standard partition of fa, b]; that is, 


Pr = {{xo, £1], (x1, x2], see tetas (3) 
where 
.b—a 
zti=za+i , fori=0,1,2,...,n. (4) 
n 


1 The Riemann integral 


Now f is increasing, so on each subinterval [x;_1, x;|, for i = 1,2,...,n, 
YA 
Mi; = f(m=i) and Mi = f(t); 
see Figure 8. Also, 62; = (b — a)/n, for i = 1,2,... ,n. Hence y= f(z) 
i=1 
b— a u = 
sm" 2_ (fai) — fle) to=@ ti T2 = b=Xp, I 
b—a = Figure 8 The standard 
TU (f(@n) — f(0)) partition for an increasing 
b-a function f 
= ZE CFO) — F(a)) 


®. Here we have used the fact that the terms of the series are 


fzr) — fleo), f(v2)— f(m); =,  f(tn)— f(tn=1), 
so everything except f(xg) and f(x,) cancels. 8 
= 


The sequence ((b — a)(f(b) — f(a))/n) is null, so it follows from 
Riemann’s Criterion that f is integrable on [a,b]. a 


Theorem F49 


A function f which is continuous on [a,b] is integrable on Ja, b]. 


Proof We use the fact that f must be uniformly continuous on [a,b]. 


®. You saw in Theorem F19 in Subsection 4.2 of Unit F1 that a function 
which is continuous on a bounded closed interval is uniformly continuous 


there. & 
Let s > 0. Then there exists 6 > 0 such that 
f(x) — f(y)| = == for all x,y € [a,b] with |x — y| < ô. (5) 


®. We use ¢/(b— a) here in order to obtain £ later in the proof. 8 


Next we choose N € N such that (b6— a)/N < ô. For n> N, let P, be the 
standard partition of [a,b] given by equations (3) and (4). 


Now f is continuous on each subinterval |1;—1, x;|, for i = 1,2,...,n. Thus, 
by the Extreme Value Theorem (see Subsection 3.3 of Unit D4), there exist 
points c; and d; in [x;-1, zi] such that 

mi = f(c) and Mi= f(d;). (6) 
Since [x;-1, xi] has length (b — a)/n < ô, we deduce by statements (5) 
and (6) that 


& 
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= 
®. Here we have used the fact that |d; — ci| < es , from which it 
n 


follows that |f(d;) — f(c:)| = Mi — mi; < — by statement (5). 8 
—a 


Hence, for n > N we have 


Thus (U(f, Pn) — L(f, Pn)) is a null sequence, so it follows from 
Riemann’s Criterion that f is integrable on [a,b]. E 


Theorems F48 and F49 show that monotonic functions and continuous 
functions are integrable, but we know that some bounded functions are not 
integrable: for example, the Dirichlet function, as you saw in Worked 
Exercise F29. This suggests the question: precisely which bounded 
functions are integrable? 


The full answer to this question is rather complicated but, roughly 
speaking, a bounded function is integrable on [a,b] if and only if it is 
continuous at ‘most’ points of [a,b]. However, it is possible for a function 
to be discontinuous at infinitely many points of [a,b] and yet be integrable 
on [a,b]. For example, the Riemann function which you met in Section 3 of 
Unit F1 is discontinuous at all rational points and yet it can be shown to 
be integrable on [0,1], the value of its integral being 0 (we do not prove 
this here). 


Riemann and Lebesgue integration 


Georg Friedrich Bernhard Riemann (1826-1866) laid down the 
fundamental ideas of the integral that is now named after him when 
he was writing his doctoral thesis in 1854. This was published 
posthumously in 1867. His formulation and proof were rather obscure 
and the version that is generally used today (including in this module) 
was given by the French mathematician Gaston Darboux (1842-1917) 
in 1875. 


In 1902 a different definition of the integral was given by another 
French mathematician, Henri Léon Lebesgue (1875-1941). The main 
difference between the two definitions is in the way the area under the 
curve is measured. The Riemann integral considers the area as being 
made up of vertical rectangles, while the Lebesgue integral considers 
horizontal rectangles. Or to put it another way, the Riemann integral 
considers the domain of the function while the Lebesgue integral 
considers the codomain of the function. Although Lebesgue’s 
definition has the advantage that it is applicable to a larger class of 
functions than Riemann’s definition, it requires the formal notion of a 
measure. (You may study measure theory in the future, but this topic 
is not covered in M208.) 


1 The Riemann integral 


Lebesgue himself provided a rather nice example to illustrate the 
difference between his approach and that of Riemann: 


I have to pay a certain sum, which I have collected in my pocket. 
I take the bills and coins out of my pocket and give them to the 
creditor in the order I find them until I have reached the total 
sum. This is the Riemann integral. But I can proceed differently. 
After I have taken out all my money I order the bills and coins 
according to identical values and then I pay the several heaps 
one after another to the creditor. This is my integral. 


(Source: Siegmund-Schultze, R. (2008) ‘Henri Lebesgue’, in Gower, T. (ed) 
The Princeton Companion to Mathematics, Princeton, Princeton University 
Press, p. 796.) 


Henri Léon Lebesgue 


1.3 Properties of integrals 


In your previous study of integration (for example, in a calculus course) 
you will have met many properties of integrals without a clear explanation 
of exactly why they hold. We now look at several of these properties and 
in some cases give an outline of how they follow from the definition of an 
integral in Subsection 1.1 and the various theorems that you met in 
Subsection 1.2. In reading this subsection, you should concentrate on 
understanding the various properties themselves; our comments on why 
the properties hold are optional reading in case you are interested. 


First, we look at the limits of the integral. We have already defined 
b 


integrals of the form f, where a < b; we now look at the situation 


a 
where a = b or a > b. 


Definitions 


Let f be a bounded function that is integrable on an interval I 
containing a and b, where a < b. Then we make the following 
definitions. 


° T= 0 
a b 
ef t=- 
b a 
0 1 
Thus, for example, f x dx is defined to equal — f x dx; you have already 
1 0 


1 0 
seen in Exercise F36 that | zdr = i, so we define the value of f x dgx to 
0 1 


1 
be —3- 
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Next we look at the integrability of a bounded function on intervals with 
endpoints a, b and c, irrespective of the order of these endpoints on the 
z-axis. Figure 9 illustrates two possibilities: a < c < b and a < b < c. The 
following result applies whatever the order of a, b and c. 


YA YA 
y = f(x) y = f(z) 


(a) (b) 


Figure 9 Integrating a function on intervals with endpoints a, b and c 
(a) when a < c < b and (b) whena<b<c 


Theorem F50 Additivity of integrals 


Let f be a bounded function that is integrable on an interval / 
containing the points a, b and c. Then 


de = A 
ey 


This result can be proved directly from the definitions above, and the 
definition of the integral given in Subsection 1.1. Notice that, in the 


b 
situation illustrated in Figure 9(b), I f is negative since b < c. 
C 


Our next result says that the integral of a non-negative function is 
non-negative, and the integral of a non-positive function is non-positive; 
see Figure 10. This result can also be proved directly from the definition of 
the integral. 
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(a) (b) 


Figure 10 The sign of the integral for (a) a non-negative function and (b) a 
non-positive function 


Theorem F51 Sign of an integral 
Let f be a bounded function that is integrable on [a,b]. 


b 
O If f(x) >0 on [a,b], then f Jf 20: 


b 
O If f(x) S U on |a|; then f Fao 


Of course, a function can be non-negative on some parts of its domain and 
non-positive on other parts. In these circumstances, subintervals where the 
function takes negative values make a negative contribution to the total 
area between the graph of the function and the x-axis. If the interval [a,b] 
on which a function is defined is made up of finitely many subintervals 
where the function is either always i a or always non-positive, 


then we can evaluate the integral J f by applying Theorem F50. If there 
a 
are infinitely many such subintervals, then the evaluation of the integral 


involves summing an infinite series; we do not pursue this here. 


Next we note that if a function f is integrable, then so is |f|; see Figure 11. 


YA YA 


(a) (b) 
Figure 11 (a) The integral of f (b) The integral of |f| 
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Theorem F52 Modulus Rule 
If f is integrable on [a,b], then |f| is also integrable on Ja, b]. 


Outline of the proof (optional) 
For a bounded function f on [a,b] and a partition of [a,b] 
PS {[xo, £1], [£1, £2], shine [Zn—-1, Ln] }, 


we define, for i = 1,2,...,n, the variation w;(f) of f over the subinterval 
[xi—1, xi] to be 


wi(f) = sup{| f(x) — f(y)|: x,y € [vi~ mil}. 
It can be shown that 

wi(f) = sup{ f(x) : x € (x;—-1,1;|} — inf{ f(a) : x € (1;—1, zil}. 
(We omit the details here.) 


Hence 
U(f,P) — L(f,P) = >) wil f)d2i, (7) 
=l 


where 6x; = £i — £i—1, as usual. 


Now, by the backwards form of the Triangle Inequality which you met in 
Subsection 3.1 of Unit D1, we have 


Uf(e)|— |f(y)|| =|f(x) — f(y)|, for x,y € (1;—1, xi], 
so that 
wi(lf|) <wi(f), fori=1,2,...,n. 
Hence, by equation (7), 
U(|fl, P)— L(F|, P) =U(f,P)— L(f, P), (8) 


for any partition P of [a,b]. We can now use Riemann’s Criterion to 
deduce from inequality (8) that if f is integrable on [a,b], then so is |f|. E 


Finally, we set out the Combination Rules for integrable functions; we use 
these a great deal to construct ‘new integrable functions from old’. 


1 The Riemann integral 


Theorem F53 Combination Rules for integrable functions 
If f and g are integrable on [a,b], then so are the following functions. 


Sum Rule f +g, with integral 


furo=fs+ fs 


Multiple Rule Af, for A € R, with integral 


Ne) 


Product Rule fg 
Quotient Rule f/g, provided that 1/g is bounded on [a,b]. 


You will meet some technigues for finding the integrals of products and 
guotients in the next section. 


The Combination Rules can be proved using the same approach as for the 
Modulus Rule. The proof uses the following inegualities, which relate the 
variations of the new functions over a subinterval (1;—1, %;| of a partition to 
those of the known integrable functions f and g: 


o wi(f +g) = wi(f) + wilg) 
e wi(Af) < |Alwi(f), for AER 
e wi(fg) = M(wi(f) + wi(g)), where M = max{sup |f], sup |g|}- 


1.4 Proofs of Theorem F44 and 
Theorem F46 (optional) 


In the proofs of these results, we will use the notion of a refinement of a 
partition. If P is a partition of an interval [a,b], then any partition 
obtained from P by adding to it a finite number of partition points is 
called a refinement of P. The partition of [a,b] obtained from two 
partitions P and P’ of [a,b] by using all their partition points is called the 
common refinement of P and P’. 
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For example, for the partitions 
P= {[0, 3l [4-3], [2.4] 15,1) 

and 
P' = { [0,5]; [5 3) 15,4} 

of [0, 1], the common refinement of P and P’ is the partition 
110,41, 4,31, [ah (2,34 1531 12113, 


as illustrated in Figure 12. 


10) En i 13 i. + 
4 2 4 
1 12 / 
0 13 13 1 
I 
y 7 y id y 
| I 
| I 
j } — / common 
0 L A 1 2 3 1 
4 3 2 3 4 refinement 


Figure 12 The common refinement of two partitions of [0, 1] 


We first prove Theorem F44, which is restated below. 


Theorem F44 


Let f be a bounded function on [a,b], and let P and P’ be partitions 


of [a,b]. Then 
UP, P) S UP): 


Proof In our proof we will assume that f is non-negative on [a,b]. The 
general result can be deduced by applying the ‘non-negative version’ of the 
result to the function g = f + c, where c is a constant so large that g is 
non-negative on [a,b]; such a constant c exists since f is bounded on Ja, b]. 


First recall that, from Theorem F43, for any partition P of [a,b] we have 


L(f,P) <U(f, P). 


We will prove Theorem F44 by showing that if P” is the common 


refinement of P and P’, then 
LP PyalG P)}SUJ, PP} SUP]: 
We claim that adding a new partition point x’ to a partition 


P= {[z0, z1], (x1, £2], sey =e En] } 
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of [a,b] does not increase the upper Riemann sum and may decrease it. 
This is because the new point 2’ lies in a subinterval [x;_1, xi], for some i, 
and the only effect on the upper Riemann sum of adding 2’ is to replace 
the rectangle with side M; standing on (x;—1,%;| with a pair of adjacent 
rectangles standing on [x;_1,x;] with heights at most M;, as shown in 
Figure 13. 


Ti—1 Ti Ti—1 x Ti 


Figure 13 The effect on the upper Riemann sum of adding a new partition 
point 


Hence adding a finite number of new partition points to P does not 
increase the upper Riemann sum. Similarly, adding a finite number of new 
partition points to P does not decrease the lower Riemann sum. 


Now the partition P” can be formed from either P or P’ by adding a finite 
number of partition points, so inequalities (9) follow. E 


To end this subsection we give the proof of Theorem F46. 


Theorem F46 


If f is an integrable function on [a,b] and (Pp) is a sequence of 
partitions of [a,b] such that ||P,|| > 0, then 


Nima URR) = lim, U Po B Í. 


n—> o0 


b 
Proof We prove that im Uf Pa) = I f. The proof that 


Jim, LAS Ps y= ff f is similar. 


We assume that f is non-negative on [a,b] so, for some M € R, we have 
0< f(x) <M, fora<2¢< bd. (10) 


Let € = 0. Since f is integrable on [a,b], there is a partition 


P= {{x0, 21), (21, 25), e.’ E E 


of [a,b] where x) = a and z}, = b, with m subintervals, such that 
b 
UEP) < f Fte (11) 


®. We use Je here in order to obtain € later in the proof. ® 
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M; 


> 


y 


7 
Tk—1 Ti Tk 


Figure 14 Removing a point 
x; of a partition 
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Now consider any partition in the sequence (P,,), of the form 
Pe = {ip Tih ertha Sol}, 

where p is the number of subintervals in Pa. For k = 1,2,...,p, we define 
My = sup{f(1): 1,1 < £ < £k} and 62, = Ek — Tk—1- 


Let P/ denote the common refinement of P’ and P,. Then, as in the proof 
of Theorem F44, we have 


UG.) < U(f,P'). (12) 
Now we can obtain P, from P’ by removing at most m — 1 of the partition 
points of P'. 


®. This is because P’ has m + 1 partition points, and all partitions have 
the points a and b in common. 8 


Removing such a point x/, lying in (xk—1, xk) say, we increase the upper 
Riemann sum by at most Mk(xx — xk—1), as illustrated in Figure 14. So, 
since M;, < M, by inequality (10), and £k — £k—1 < ||P,||, for 

k=l 2.20250; 


U(f, Pa) < U(f, Pa) + (m—1)M||Pall. (13) 
Combining inequalities (11), (12) and (13), we obtain 


b 
(SPa) < (| f+ 3e) + (m-DMIPaI. 
a 
Since ||P,,|| + 0 as n + co, we can choose N so large that 
(m—1)M||Pal| < $£, for alln >N. 


®. We use Zé here in order to obtain € in the next step of the proof. ® 


Hence 


b b 
U(f,Pi) £ (/ f+4e] +je= | fte, foralln>N. 
b 
Since U(f, Pa) 2 f f, by the definition of the integral, we deduce that 


b b 
vue- f f =U P) - f aa hainesi 


a 


b 
®. So U(f, Pn) — f f tends to zero as n tends to co. .& 


b 
Hence lim U(f, Pa) = / f, as required. Oo 
n—> o0 a 


2 Evaluation of integrals 


In Section 1 you saw what it means for a bounded function defined on a 
closed interval to be integrable. In this section you will study a variety of 
techniques for evaluating the integrals of such functions. 


2.1 The Fundamental Theorem of Calculus 


In this subsection we show that integration and differentiation are 
intimately related by proving a result known as the Fundamental Theorem 
of Calculus. You probably know this result from a previous course on 
calculus, but it is worth pausing to reflect on how remarkable it is. In 
Section 1 we defined the integral by means of more and more accurate 
estimates of the area between the graph of a function and the x-axis. At 
no point did it seem that this process was related to differentiation, but 
the Fundamental Theorem of Calculus shows that integration and 
differentiation are in some sense inverse processes. As you will see, this 
fact is enormously helpful in evaluating integrals. 


We begin our exploration of these ideas by defining a primitive of a 
function. 


Definition 
Let f be a function defined on an interval J. Then a function F is a 
primitive of f on I if F is differentiable on J and 


PQ) = ja), termes. 


It follows from this definition that finding a primitive is the inverse of 
finding a derivative. (For this reason, a primitive is sometimes called an 
antiderivative.) Note that the domain of the primitive F may be larger 
than the interval J on which f is defined. 


As an example of finding a primitive, let 
f(x) = tanz. 

Then the function 
F(x) = log(sec x) 


is a primitive of f on the interval (—7r/2, 7/2), since 


1 
F"(x) = = secrtane = tang, for x € (—a/2,7/2). 
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Exercise F38 
(a) Let 
f(a) = (z? — 4) 72 (w € (2, 00)). 
Prove that 
F(a) = log (z + (z? — 4)/?) 


is a primitive of f on (2,00). 


(b) Let 
f(z) = sech (= È ) 


cosh x 


Prove that 
F(x) = tan7*(sinh z) 


is a primitive of f on R. 


We now state and prove our main result, the Fundamental Theorem of 
Calculus. The theorem tells us that we can evaluate the integral of a 
function f on an interval [a,b] by finding a primitive F of f on [a,b]. Note 
that the expression F(b) — F(a) is sometimes written as (F(x)|? Or F(2)|°. 
Theorem F54 Fundamental Theorem of Calculus 


Let f be integrable on [a,b] and let F be a primitive of f on [a,b]. 
Then 


b 
f f=FO-FO). 


Proof Let 

Pa = {(x0,2%1),.+-, |2;—1,%;),+++,|2a—1,2n|}, forn =1,2,..., 
be a sequence of partitions of [a,b], with zo = a, £n = b and ||P,|| > 0. 
On each subinterval 

12;=1,0;|; fori =i; eeii 


the function F satisfies the conditions of the Mean Value Theorem which 
you met in Subsection 4.1 of Unit F2 Differentiation, since a primitive is 
differentiable and hence continuous. Thus there exists a point 
Ci E€ (£i—1, zi) such that 
F (ai) — F(xi-1) = F' (ci) (xi — 271) 
= f(ci) xi, (14) 


where 62; = £i — £i—1. 


@. Now recall that we use the notation m; = inf{f(x) : x € (1;—1, %;|} and 
M; = sup{ f(x): z € [ai-1,0;]}. © 


Since 
mis Te) My tora =— 1,2, 410.5% 


it follows that 
i=1 i=1 i=1 


Using equation (14), we can rewrite this statement as 
n 
Les) a a aU). 
i=1 


®. We now use telescopic cancellation to evaluate the series in these 
inequalities, which is 

(F(a1) — F(zo)) + (F(22) — F(21)) +++ + (F(tn) - F(@n-1)). 
The only remaining terms after cancellation are F'(x,)) and —F (zo). ® 
The series has sum F'(x,) — F(xo) = F(b) — F(a), so 

L(f, Pn) < F(b) — F(a) < UCF, Pn). (15) 


Since f is integrable on [a,b], the sequences (L(f, P,)) and (U(f, Pn)) both 
b 
converge to f f, by Theorem F46. It follows from inequalities (15) and 


the Limit Inequality Rule for sequences (Theorem D11 in Unit D2 
Sequences) that 


[rsr@-re@s fs, 


which gives the required result. E 


b 
Theorem F54 shows the close relationship between the integral f f and 
a 


any primitive F of f on [a,b]. Because of this result, a primitive F is also 
called an indefinite integral of f and denoted by f f(x)dx. Moreover, 


the process of finding a primitive of f is often informally called 
integrating f, and in this context the function f is called an integrand. 


b 
Also, the integral f f is often referred to as the definite integral of f 
over [a,b]. g 


We can use Theorem F54 and the table of standard primitives at the end 
of this unit to evaluate many integrals. 
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Worked Exercise F30 


1 
Evaluate J 2” dr. 
0 


Exercise F39 


Using the Fundamental Theorem of Calculus and the table of standard 
primitives, evaluate the following integrals. 


(a) [Eta (b) f oszd 


2.2 Primitives 


It is natural to ask: can a function have more than one primitive on an 
interval? The answer to this question is ‘yes’: for example, on (—1, 1) the 
functions 


rm r? and ze r? +1 
are both primitives of the function 
GH 27. 


However, any two primitives of a function f on an interval can differ only 
by a constant. 


Theorem F55 Uniqueness Theorem for Primitives 


Let Fy and F be primitives of f on an interval J. Then there exists 
some constant c such that 


F(x) = Fi(z) +c, forge. 


Proof Since F} and F, are primitives of f on J, 
Fi(z)= f(x) and F}(x)= f(x), forz el, 
SO 
F(x) — F{(x) =0, forz eI. 


Thus, by the Zero Derivative Theorem (Corollary F39 in Unit F2), there 
exists a constant c such that 


F(x) _ F(x) =c, forrE€ 1. a 


The range of primitives we can find is considerably extended by the use of 
several Combination Rules. T'hese rules can be proved using the 
corresponding rules for derivatives; we omit the details. 


Theorem F56 Combination Rules for primitives 


Let F and G be primitives of f and g, respectively, on an interval J, 
and let A € R. Then, on I: 


Sum Rule f +g has a primitive F + G 
Multiple Rule Af has a primitive AF 
Scaling Rule x> f(z) has a primitive x +> Z F(a), 
tone UE 


For example, it follows from the table of standard primitives and the 
Combination Rules that the function with domain R* and rule 


at 3x7! + sinh 2x 
has a primitive 
x e 3log x + 5 cosh 2x. 


In applications of these Combination Rules we do not usually mention the 
rules explicitly. 


Exercise F40 


Using the table of standard primitives and the Combination Rules, find a 
primitive of each of the following functions. 


(a) f(x) =4logz —2/(4+ 27) (x € (0,0)) 
(b) f(z) =2tan3z+e"cosx (x € (—1/6,7/6)) 
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2.3 Techniques of integration 


The Fundamental Theorem of Calculus provides a powerful method for 
evaluating certain integrals. However, even when we know that a function 
f has a primitive F, it may not be possible to determine F explicitly. 


In fact, most functions f, even quite simple ones, have primitives which are 
not standard functions. For example, the primitives 


s dx 


cannot be expressed as a combination of a finite number of rational, nth 
root, trigonometric, exponential and logarithmic functions. (Note that 


dx 
here we have used a standard shorthand and written i ——— instead of 


I (log x)? 
J (loga)? 12) 


However, by using the Combination Rules we can certainly integrate any 
polynomial function, and the primitive is then always another polynomial 
function. For example, 


(2? 045) de = 15 — $a? + 50. 


There is also is a standard procedure called the method of partial fractions 
for integrating a rational function. This procedure is not used in this 
module but is often used for evaluating the integrals of rational functions in 
complex analysis. Using this method, it can be shown that any primitive of 
a rational function can always be expressed in terms of rational functions, 
logarithms of rational functions and inverse tangents of linear functions. 


There are also various techniques which can be applied to certain other 
types of function; the art of integration lies in recognising these types. We 
now describe briefly some basic techniques of integration. 


Integration by substitution 


We describe two related techniques of integration by substitution. The 
first is used when the integrand is of the form 


zr— f(9(x))g'(2). 
In this case, if F is a primitive of f, then 


d 


Tat a2) = F(g(r))g'(z) = F(g(z))9"(@), 


by the Composition Rule for derivatives (also called the Chain Rule). Thus 
x > F(g(x)) is a primitive of  —> f(g(x))g’(x) on any interval in the 
domain of F o g. So if we substitute u = g(x), then we obtain a simpler 
integrand since 


ECOG dr = F(g(x)) = F(u) = fro du. (16) 


This technique is worth trying if you can express the integrand in the form 
f(g(a))g'(x), for some functions f and g, as illustrated in the next worked 
exercise. 


Worked Exercise F31 


Find a primitive of the function 


zr — (2 +1)? (TER). 


Solution 
®. Differentiating x? + 1 gives 3x7, which suggests that we try putting 
u=g(z)=2°+1. SE 
d 
Put u = 9(z) = 2° +1. Then = = gf e) = 30,60 CM = SV" da 
Ji 


®. Now substitute to express the integral in terms of u, adjusting by 
a multiplicative constant as necessary. 8 


We have 


NG + de = if + 1)° 327 dx 


SE du. 


@. Next, evaluate the integral in terms of u. 8 


EGO 
x gu 


®. Finally, substitute back for u to obtain the result in terms of x. ® 
= 4(2° +1)’. 


Thus +(x° + 1)° is a primitive of a (e + 1)8. 
27 


®. You can always check your result by differentiating. © 


The following strategy summarises the approach taken in Worked 
Exercise F31. 
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Strategy F8 


To find a primitive | f(g(x))g'(x) dx using integration by 


substitution, do the following. 


d 
1. Choose u = g(x). Find = = g'(x) and hence express du in terms 
i 


of x and dz. 
2. Substitute u = g(x) and replace g'(x) dx by du (adjusting constants 


if necessary) to give | f(u) du. 


3. Find fro du. 


4. Substitute u = g(x) to give the required primitive in terms of x. 


Remarks 


1. If we are evaluating an integral, rather than finding a primitive, then 


there is no need to perform step 4 of Strategy F8. Instead, we can 
change the x-limits of integration into the corresponding u-limits: 


b g(b) 
f CEOE / Hida 
a g(a) 


For example, in Worked Exercise F31 we have 
u= g(x) = 2° +1, 
SO 
whenx=0, w=1, andwhenz=1, u=2. 


Hence 


1 2 
) r? (x3 + 1)8 dx = t f u® du. 
0 1 


. When applying equation (16), you may be able to spot a primitive F 


of f immediately, in which case you can write down the required 
primitive F(g(x)) directly without need for a substitution. For example, 
if you were able to spot the primitive in Worked Exercise F31, you 
could just write 


foe +1) dz = fe + 1)? 8a? dx 


l 
3 
ix $ (2° +1)? = (2 +1)°. 


One particular situation in which this simplification can often be 
applied is with an integrand of the form g'(x)/g(x) on an interval J, 
since we have 


[2 dz = log(g(x)), if g(x) = 0, forzre J, (17) 
g(x) 
because < log(g(x)) = aa x g'(x): 


An example involving this useful formula appears in the next exercise. 


Exercise F41 


Find a primitive of each of the following functions. 
(a) f(x) =sin(sin3z)cos3x (a €R) 

(b) f(x) =27(2+323)" (TER) 

(c) f(z) =asin(2x”) («€R) 

(da) f(z) =2/(2+327) (TER) 


T 


Exercise F42 


Evaluate the integral 
1 x 
e 
— d. 
J +e” 


Our second substitution technique is a modification of the above method, 
which we call backwards substitution. It is based on the formula 


reg de = f EDK du, 


obtained from equation (16) by swapping the variables x and u. Here h is 
a function such that x = h(u), often found by first writing u = g(x) 
where g has an inverse function g~! = h. Backwards substitution is worth 
trying if it makes part of the integrand significantly simpler. The next 
worked exercise gives an example of the use of this technique. 
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Worked Exercise F32 


Find a primitive of the function 


e2t 


opie CE (0.00)). 


ut 


@®. If u = g(x) = (e? — 1)'/2, then g has an inverse function 
x = h(u) = log(u? + 1). # 


Put u = g(x) = (e? — 1)!/2. Then the inverse function is 


d 
= U = log(u? + 1), so = = pu) = 7 = i and therefore 
u 
®. Now substitute to express the integral in terms of u. ® 
We have 


e27 elog(u? +1)? 2u 
sihr =} — d 
Ne ie) 


EBI) 2 
- [H = S 
u u? + 1 


= fre + 1)du 


@. Next, evaluate the integral in terms of u. © 


= zu) ap ZU 


$. Finally, substitute back for u to obtain the result in terms of x. 8 


= le" — 1)9/? + (e — 1)”. 


e2t 


Tus (e" = IVE ME" = EA a primitive at, Cae 
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The following strategy summarises the approach taken in Worked 
Exercise F32. 


Strategy F9 


To find a primitive J f(x) dx using integration by backwards 
substitution, do the following. 
1. Choose u = g(x), where g has an inverse function x = h(u). Find 


as 
a h'(u) and hence express dx in terms of u and du. 
u 


2. Substitute x = h(u) and replace dx by h’(u) du to give a primitive 
in terms of u. 


3. Find this primitive. 


4. Substitute u = g(x) to give the required primitive in terms of z. 


As before, if we are evaluating an integral, then instead of step 4 in 
Strategy F9, we can change the z-limits of integration into the 
corresponding u-limits: 


b (b) 
f ieas i ~ O a 
a g(a) 


For example, in Worked Exercise F32 we have 


u=g(z) = (e*- 1)", 


SO 
when z = log2, u = 1, and when zx =log3, u = V2. 
Hence 
log 3 e27 v2 2 
mt] 2(u + 1) du. 
J (e7 = 1)1/2 1 


Exercise F43 
(a) Find a primitive of the function 


f(x) = : 


32-13 +2(2—)? (x € (1,00)), 


using the substitution u = (x — 1) 72. 


(b) Evaluate the integral 


log 3 
f e” y1 + e dz. 
0 
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Integration by parts 


The technique of integration by parts is derived from the Product Rule for 
differentiation, 


(fo = F'g + Fg', 
which implies that 


fid = to- | Pa so [setae - f ta 


This formula converts the problem of finding a primitive of fg’ into the 
problem of finding a primitive of f'g. Integration by parts is worth trying 
if you can express the integrand as a product of two functions, f(x)g'(x), 
where f(x) becomes simpler on differentiation, and g'(x) becomes not 
much more complicated on integration. 


Here is a strategy for using integration by parts. 


Strategy F10 


To find a primitive | k(x) dx using integration by parts, do the 
following. 


1. Write the original function k in the form fg’, where f is a function 
that you can differentiate and g’ is a function that you can 
integrate. 


2. Use the formula fr = fo- | ts. 


The next two worked exercises give examples of using this technique. 


Worked Exercise F33 


Find 


fe cos x dx. 


Sometimes we have to multiply the integrand by the factor 1 in order to be 
able to apply integration by parts, as in the following example. 


Worked Exercise F34 


Evaluate the integral 


1 
J tan! zdz. 
0 


Exercise F44 


(a) Find a primitive of the function 
k(x) =a loge (a € Rt). 
(b) Evaluate the integral 


T/2 
2 
J xi cos z dx. 
0 


Hint: Use integration by parts twice. 
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Sometimes we need to evaluate an integral [,, that involves a non-negative 
integer n. A common approach to such integrals is to relate the value of In 
to the value of In—1 or In_2 by a reduction formula (sometimes called a 
recurrence relation) using integration by parts. Here is an example that 
will be important later in the unit. 


Worked Exercise F35 
Let 
a /2 
n=} sin” gdr, n=0,1,2,.... 
0 


(a) Evaluate Jp and 14. 
(b) Prove that 


Ma = (" — *) In-2, forn> 2. 
n 
(c) Deduce the values of I>, 13, I4 and Is. 
Solution 
an /2 
(a) We have w= Lan = m2 emal 
0 
T/2 = 
u) Ssmsd = [= csn] = 

0 

(b) We write 


a /2 
Iy = U sin x sin”! g dz. 
0 


@. We integrate sin x and differentiate sin" 12. © 


Using integration by parts we find that, for n > 2, 


In = [(— cos 2) sin"! a] ue = 


a /2 
| (— cos x)(n — 1) sin" 2? z cos x dz 
0 
a /2 
=0+(n-1) | cos? z sin”? z dx 
0 


a /2 
= (n— 1) | (1 — sin? 1) sin" ? zdr 
0 


T2 T2 
= (n— 1) sin”? z dx -f sin” x dx 
0 0 


= (n—1)Un-2 — In). 
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3 Inequalities, sequences and series 


By repeatedly applying the reduction formula in Worked Exercise F35, we 
obtain the general formulas 


| ae eae 

MDTA G In ě 2 
and 

E. 5. A eee _2n 

enpi g S T7 mL 


We will use these formulas in Subsection 3.2. 


Exercise F45 
Let 
1 
de =| CL dr, awa eek 
0 
(a) Evaluate Jp. 
(b) Prove that 
In =e-—nIn1, forn=1,2,.... 


(c) Deduce the values of I, I2, 13 and 14. 


3 Inequalities, sequences and series 


Often it is not possible to evaluate an integral explicitly, and a numerical 
estimate for its value is sufficient. This situation can arise both in 
applications of mathematics and in proofs that involve integration. In this 
section we study some inequalities satisfied by integrals, and apply these to 
find two remarkable formulas for 7, and to decide whether certain series 
are convergent or divergent. 
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3.1 Inequalities for integrals 


The basic inequality rules for integrals are as follows. 


Theorem F57 Inequality Rules 
Let f and g be integrable on Ja, b]. 
(a) If f(x) < g(x), for x € [a,b], then 


b b 
oa 
(b) Ifm < f(z) <M, for x € [a,b], then 
b 
mo-a) < | o ey. 


Proof The proofs of parts (a) and (b) of Theorem F57 are illustrated in 
Figure 15. 


YA 


(a) 


Figure 15 The proofs of the Inequality Rules parts (a) and (b) 
(a) Let P be any partition of [a,b]. Since 
f(x) < g(x), for x € [a,b], 


the infimum of f on each subinterval of P is less than or equal to the 
infimum of g on that subinterval, and therefore L(f, P) < L(g, P). 
Thus 


b b 
/ fail. P)< mp l@,P)= | " 
a P P a 


since f and g are both integrable on [a, b|. 


(b) Since f(x) < M for x € [a, 6], it follows from part (a), with g(x) = M, 
that 


[ors f mar= Moa 


The proof of the left-hand inequality is similar. | 
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The Inequality Rules allow us to estimate a complicated integral by 
evaluating a simpler one, as in the next worked exercise. 


Worked Exercise F36 


Prove the following inequalities. 


1 x 3 
——___ dr< ! log2 b =<] — 
(a) I 2— sinf r "SAE (b) v34 J-1V¥24+ 2° 


2 
dx <3 


Solution 


(a) 


Since 

|sin z| < |x|, forx ER, 
®. This is Corollary D46 from Subsection 2.3 of Unit D4. & 
it follows that 


sinf x < r$, for x eR. 


Thus 2 ~ sin q > 2-—2*>0 for x € [0,1], so 


z? 


— — fi Ome 
S=sn 2 2— 7°" ee Ue 


Hence, by Inequality Rule (a), we have 


1 3 1 3 
dg x 
o 2 = ginal ge 0 2=e 


1 
= [—4 log(2 — ae 
= —}(log 1 — log 2) 
= ; log De 


@. Here we have used equation (17) with g(x) = 2 — zt. ® 
Since the function x —> v2 + 2° is increasing on [—1, 2], we have 


1 < V¥2+2° < V34, for x € [-1,2], 
so 
1 il 
= 3 SS 
V34 v2+r 
Since the length of the interval [—1, 2] is 3, it follows from 
Inequality Rule (b) that 


L does e [=LA] 


2 
—<| mm 
V 34 -1V2+2° 
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Exercise F46 
Use the Inequality Rules to prove the following inequalities. 


3 1/2 
a xsin(1/x1?) dx < 4 b) iz e? dx < le!/4 
1 ? 0 s 


We saw in Subsection 1.3 that if the function f is integrable on fa, b], then 
so is |f|. We now use the Inequality Rules to obtain an inequality 
involving the integrals of f and |f|, known as the Triangle Inequality for 
integrals. The name arises because of the similarity between this inequality 
and the Triangle Inequality for numbers: 


n n 
2 a| < $ lail. 
i=1 i=1 


Theorem F58 Triangle Inequality for integrals 
Let f be integrable on [a,b]. Then 


iG < fit. 


Furthermore, if |f (x)| < M for x € [a,b], then 
b 
ie 
Proof We know that, for all z € [a, b|, 


-|f| < f(x) < |f). 
Since |f| is integrable on [a,b], it follows from Inequality Rule (a) that 


- fis firs fa 


which is equivalent to 


fi < fish 


Finally, if |f (x)| < M for x € [a,b], then, by the above inequality and 
Inequality Rule (b), 


E 


<M(b-a). 


b 
< | |f| < M- a). 


3 


Sometimes we use the Inequality Rules and the Triangle Inequality in 
combination, as in the next worked exercise. 


Worked Exercise F37 


Prove that 


Solution 


By the Triangle Inequality for integrals, 


T/2 = T2 = 
i x—m/2 T < [ x— 7/2 
0 2+cosz 0 22 COS GE 
T/2 Ga 
=| T/2-2 ix 
0 24,cosx 


@. For 0 < x < 1/2, we have r/2— x > 0 and cosx>0. # 


Next, since 
2 +cosr2 2, for x € [0,7/2], 


we have 
1 


1 
ez a 
Sen or eee 


Thus, by Inequality Rule (a) and statement (+), 


ee Na 
| Diao i; 3(5- 2) a 


IA 


U 
3x NIK DI a 


je 
o>) 


Exercise F47 


Prove the following inequalities. 


4 sin(1/2) 
fi) | =t (b) 


v tanx 
——— dx 
o 3-—sin(x?) 
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3.2 Wallis’ Formula 


In Worked Exercise F35 we used a reduction formula to show that if 


a /2 
as, sid” tdr, 9 SN 1 2s. 
0 


then 
—1 

Mimi Ij=1 and [,= (" ) Inn for n 2 2. (18) 

We also remarked that by repeatedly applying eguations (18) we obtain 
1 3 5 2n—1 7 

Jon =] 2424 Be ee Se, OF Si 19 

an 2°46 i U oe uu 
and 

246 2n 
Í ESUE E Ff > 1. 20 
mam T al NE (20) 


Note that the formula for Jz, involves 7, but the formula for [9,41 does 
not. 


We now use these results, together with inequalities between various 
integrals of the form In, to establish two remarkable formulas for 7, the 
first of which is Wallis’ Formula. 


Theorem F59 


; 2 2 A ad G 2n 2n T 

(a) lim at a e ee OS aa a a E E = — 

nm>o\l 3 3 5 5 7 2n—1 2n+1 2 
(Ge 2e 


Oa a 


John Wallis (1616-1703) was the most influential English 
mathematician before the rise of Isaac Newton. He was a 
cryptographer to the Parliamentarians during the Civil War and in 
1649 was appointed Savilian Professor of Geometry at Oxford. 


His most important work is his Arithmetica Infinitorum (Arithmetic 
of Infinitesimals), published in 1656. It was in this work that Wallis 
derived the formula which bears his name, and it was through 
studying this work that Newton came to discover his version of the 
binomial theorem. 


Wallis’ treatise on conic sections, published in 1655, contains the first 
publication of the symbol for infinity, oo. 


Part (c) of the next exercise gives an identity needed in the proof of 
Theorem F59. If you are short of time, then you may wish to skip this part 
of the exercise, which is quite challenging, and skim read the solution and 
also the proof of the theorem. 


3 Inequalities, sequences and series 


Exercise F48 


Forn = 1,2,..., let 
2244 6 6 2n 2n (nl)? 22” 
me eee ee + eee and bp = A, 
133 5 5 7 2n—1 Intl (Qn)l/n 
(a) Evaluate a, and bn, for n = 1, 2,3. 
(b) Verify that 


2 1 
p= — ) an, for n = 1, 2,3. 


an 


(c) Prove that 


2 1 
= ( us ) an form = 152,522. 


n n 


Proof of Theorem F59 
Let 


T/2 
p= sin” gdz, n=0,1,2,..., 
0 


and let the sequences (an) and (bn) be as given in Exercise F48. 


(a) Using equations (19) and (20), we obtain, for n > 1, 


Tmo 1-3+3-5- +++ -(Qn-1W(Qn+1) er 1 2 
Iny 2 2+ 4+ 4+ +++ + (2) (2) De a 
SO 
a Iny \ T 
" Ion, J 2 
Thus to prove part (a), it is sufficient to show that 
I 
ant? +1 asn— oo. (21) 
Tan 


We do this as follows. Since 0 < sina < 1, for x € [0,7/2], we have 


2 2n+1 


sin” x > sin z 2 sin?"2 r, for x € [0,7/2]. 


It follows by Inequality Rule (a) that 
Ton 2 I2n41 2 Lanse. 
Thus, on dividing by Jon, we obtain 


1> Tanti = Ton+2 _ 2n+1 _ ey 
= lo, ~ dan 2n+2 2+2/n 


by equation (18). On taking the limit as n — co, we deduce that 
statement (21) holds, by the Squeeze Rule for sequences from Unit D2. 
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(b) We know, from Exercise F48(c), that 


Qn +1 1 
2 = ( eas ) ay = (242) an 
n n 


By part (a), 


T 
an > 7 asn > 00, 
so, by the Product Rule for sequences, 
T 
bn 2X ET as n —> OO. 


Hence, by the continuity of the square root function, 


bn > VT as n > oo. u 


3.3 The Integral Test 


In this subsection we introduce a method based on integration for 
determining the convergence or divergence of certain series of the 
0O 
form ) : f(n), where the function f is positive and decreasing and 
n= 


tends to 0. The method is based on the fact that it is often easier to 
evaluate an integral than a sum which has a similar behaviour. 


Theorem F60 Integral Test 


Let the function f be positive and decreasing on [1,0o), and suppose 
that f(z) > 0 as x > co. Then 


f(n) converges if the sequence ( | U is bounded above 
1 


f(n) diverges it [ f > œ as n > oo. 
1 


Remarks 
1. The Integral Test is also called the Maclaurin Integral Test. 


2. In both parts (a) and (b), the number 1 can be replaced by any positive 
integer. 


3 Inequalities, sequences and series 


Colin Maclaurin (1698-1746) was a Scottish mathematician who spent 

most of his career as professor of mathematics at Edinburgh 

University, having been appointed on the recommendation of Isaac 

Newton. A popular teacher, he was described as a ‘favourite professor’ 

and the ‘life and soul’ of the university. He is notable for having | 
extended Newton's work on the calculus and geometry. His 
two-volume Treatise of Fluxions (1742), which was the first systematic 
treatment of Newton’s methods, contains a detailed discussion of 
infinite series and includes the Integral Test in verbal form. 


The Integral Test was rediscovered by Cauchy — he published it in 
1827 — and consequently it is also known as the Maclaurin—Cauchy 
Integral Test. 


Colin Maclaurin 
Proof of Theorem F60 For n=2,3,..., let 
Sn = f(1) + f(2) +---+ f(n) be the nth partial sum of the series 5 f(n) 


n=1 
and let P,-1 be the standard partition of [1, n] with n — 1 subintervals of 
length 1, that is, 


{[1,2],...,/é,¢+1],...,[n-—1,n]}. 

Since f is decreasing on [1,00), we have, for i = 1,2,...,n — 1, 
mi = f(i+1) and Mi= f(i). 

This is illustrated in Figure 16. 


Ya YA 


area = M; 62; 


= f(t) 


&Y 
p 
N 
S 
>. 
+ 
= 
3 

SY 


4141 en 
Figure 16 ‘The lower and upper Riemann sums of f 


Also, each subinterval in the partition has length 1. Hence the lower and 
upper Riemann sums for f on [1,n] are 


L(f,P, n—1) -Emx 2)+}x" -+ f(n) = sn- f(1) 
and 


U(F Pai) -Fma 1) +--+ f(n- 1) = sn — f(n). 
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Since f is bounded and monotonic on [1, n], it follows from Theorem F48 


n 
in Subsection 1.2 that the integral [, = f f exists and satisfies 
1 


L(f,Pn—-1) < In < U(f, Pat), 
so 

5n — f(1) < In < sn — f(n). (22) 
We now consider two cases separately. 
Case (a): The sequence of integrals (In) is bounded above. 


In this case, there exists M € R such that 
In <M, for n= 2,3,.... 

It then follows from the left-hand inequality in statement (22) that 
sa = JU) +My tor R= 2, 3 


Thus the increasing sequence (s,) is bounded above, so it is convergent, by 
the Monotone Convergence Theorem (Theorem D22 in Subsection 5.1 of 
Unit D2). 


x oo . 
Hence the series ) ; f(n) is convergent. 
n= 


Case (b): The sequence (In) is not bounded above. 


The sequence (In) is increasing, since 


n+l 
Ini —a= fo $20, 
n 
so in this case 
In 2 œ asn 3 oo. 
It follows from the right-hand inequality in statement (22) that 
Sy. Ins for N= 25352644 


Thus, by the Squeeze Rule for sequences which tend to infinity 
(Theorem D18 in Subsection 4.3 of Unit D2), 


Sn > œ as n — OO. 


Hence the series ane i f(n) is divergent. E 
n= 


In Theorem D33 of Unit D3 Series you saw that the basic series 
IM {! /n? converges for p > 2 and diverges for 0 < p < 1. We can now 
n= 


use the Integral Test to deduce the behaviour of this series for all p > 0. 
We first consider the case when p = 1. 


3 Inequalities, sequences and series 


Worked Exercise F38 


dx 
Use the fact that J= = log x to prove that D is divergent. 
z 


n= r 


Solution 
Let 


f(a)=—= (e€ [1,00)), 
Then f is positive and decreasing on {1, co), and 
f(x) 30 asx ow. 


Also, for n > 1, 


fi he 


[ log s)" 


= logn > œ asn— oo. 


Hence, by part (b) of the Integral Test, the series diverges. 


We now consider the series i /n? when p > 0 but p #1. 
n= 


Worked Exercise F39 


Use the Integral Test to determine the behaviour of the series 


OO 


1 
ie. forp>0, pA. 


n=l 


Solution 
Let p > 0 and p £ 1, and let 


JE) = Wak (eea 
Then f is positive and decreasing on |1, co), and 
f(x) 30 as z> o. 
Also, for n € N, 
n n i= 
Sea 


First suppose that p > 1. Then p— 1 > 0, so equation (+) gives 


jz E pees PE 
1“ p= NT) ~ p-1 


| -r 
l-p 
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Exercise F49 


Show that 


f = S = 
z(logx)? log’ 
and hence prove that 


co 


1 
`> E is convergent. 


n=2 


4 Stirling’s Formula 


This section concerns the value of the quantity n! which arises in many 
problems in probability. You will see that integration techniques give an 
excellent estimate for n!, called Stirling’s Formula, which can be expressed 
as 

n! 


lim —————— = 1 
noo V2rn (n/e)” 


4.1 Comparing functions of n 


In this subsection we define a relation ~ on the set of all positive functions 
with domain N to provide a way of comparing the behaviour of such 
functions for large values of n. This relation is, in fact, an equivalence 
relation, although we do not prove this here or make explicit use of these 
properties. (You met equivalence relations in Unit A3 Mathematical 
language and proof and you may like to try to prove that ~ is an 
equivalence relation yourself; it is included as an exercise in the additional 
exercise booklet for this unit.) 


Definition 
For positive functions f and g with domain N, we write 
f(n) = g(n) as n 7 co 


to mean 


f(n) 


—— 31 anon. 


g(n) 


For example, 


n2 tnan? as n + 00, 


since n? +n > 0 and n? > 0, for n = 1,2,..., and 
n? +n 
ne 
Note that the statement 


—1 asn 0. 


f(n) = g(n) as n —> co 


does not imply that f(n) — g(n) tends to zero or is even bounded. For 
instance, in the above example, (n? + n) — n? = n tends to infinity. Note 
that often we omit ‘as n — oo’ when writing expressions of this type. 


We have the following Combination Rules for ~. These rules follow from 
the Combination Rules for sequences; we omit the proofs. 


Theorem F61 Combination Rules for ~ 
If fi(n) ~ gi(n) and fa(n) ~ ga(n), then: 

Sum Rule film) + faln) = gi(n) + ga(n) 
Multiple Rule Afi(n) = Agi(n), for A € Rt 
Product Rule fi (n) fo(n) = gi(n)ga(n) 


fi(n) _ aln) 


Quotient Rule AT 
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Note that if l > 0, then the statements 
f(n) 391 as n oo 

and 
f(n)~l asn — oo 


are equivalent, since each is equivalent to the statement 
f(m) 


= l as n — oo. 


4.2 Calculating factorials 


For small values of n we can evaluate n! directly by multiplication or by 
using a scientific calculator. 


Exercise F50 


Complete the following table of values of n!. 


n| on! n n! n n! 

1 6 720 || 20 | 2.432... x 1018 
2 7 5040 || 30 

3 8 40320 || 40 

4| 24| 9| 362880 || 50 

5 


120 || 10 | 3628800 || 60 


As n increases, n! grows very quickly, and is soon beyond the range of a 
calculator. Many calculations in probability theory involve n! for large 
values of n, so it is important to be able to estimate this quantity as 
accurately as possible. The following result gives us a way of doing this. 


Theorem F62 Stirling’s Formula 
n! ~ V2an(n/e)” asn— oo. 


We give the proof of this result in Subsection 4.3. 


James Stirling (1692-1770) was a Scottish mathematician who made 
many contributions to analysis. He was friends with Isaac Newton 
and Colin Maclaurin, and at Maclaurin’s request he provided help and 
criticisms when Maclaurin’s Treatise of Fluxions was in proof. 


Stirling’s most important work was his Methodus Differentialis, 
published in 1730. It focuses largely on infinite series and was directly 
stimulated by earlier work of Wallis and Newton. It includes a series 
for log n!, now usually written in the form known as Stirling’s 
Formula, given above. Earlier in 1730, Abraham de Moivre 
(1667-1754), famous for his work on the laws of chance, had published 
a similar formula for log n!, but he had not been able to find a precise 
value for the constant which Stirling showed to be V2r. 


Exercise F51 


Use a calculator to evaluate V2rn (n/e)" for the following values of n. 
(a) n=5 (b) n=10 (c) n= 50 


As you saw in Exercise F51, even for small values of n, Stirling’s Formula 
gives reasonable approximations to n!, and the relative error (that is, the 
error expressed as a percentage of the true value) decreases as n increases. 


n n! Stirling’s approximation Relative error 

10 3628800 3 598 696 0.83% 

20 2.433 x 1018 2.423 x 1018 0.42% 

50 3.041 x 10% 3.036 x 1064 0.16% 
100 9.333 x 10457 9.325 x 10157 0.08% 


In fact it can be shown by a more careful argument that 


el/(12n+1) < ua < e/l2n) for n > 1. 


V2rn (n/e)” 
For example, if n = 10, then 


el/(12n+1) — 91/121 — 199829... and e!/(12") = ¢1/120 _ 1.00836... 


which indicates a relative error of about 0.8%, as shown in the above table. 
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Stirling’s Formula can be used to give estimates of probabilities. For 
example, if we toss n fair coins, then it can be shown that the probability 


1 ! 
of obtaining exactly r heads is n —, where my). __™ is the 
ry 2" r ri(n—r)! 


binomial coefficient which you met in Subsection 3.4 of Unit D1 (we do not 
prove this here). Therefore, if we toss 200 coins, then the probability of 
obtaining exactly 100 heads and 100 tails is 


Ge 1 200! 
100) 2200 ~ (1001)2 2200 

V 4007 (200/e)?°° 
(V200% (100/e)190)* 9200 
4007 (20020 /e200) 


~ 3007 (100200 /¢200) 2200 


_ V4007 


2007 
1 
10./7 

1 


— 17.724...’ 
perhaps rather higher than you might expect. 


~~ 


The following exercises give you a chance to practise using Stirling’s 
Formula. In some of these exercises you will need to use the Combination 
Rules for ~ from Subsection 4.1. 


Exercise F52 


Use Stirling’s Formula to estimate each of the following numbers (giving 
your answers to two significant figures). 


300\ 1 300! 1 
a WE (b) (1001)3 3300 


Exercise F53 


Use Stirling’s Formula to determine a number A such that 


& J r) ~ 22” as n— o. 
2n n 


Exercise F54 


Use Stirling’s Formula to prove that 


Hint: You can assume that if f(n) ~ g(n), then (f(n))!/" = (g(n))!/”. 
(This result about ~ holds because if n € N, then 


1/n 


a<a/" <1, for0O<a<1l, 


and 
1<a” <a, fora>l, 


by the rules for inequalities.) 


4.3 Proof of Stirling’s Formula (optional) 


The proof of Stirling’s Formula is quite long. If you don’t have time to 
study it properly, then you may find it interesting to skim through this 
subsection and see how integration is used to obtain an estimate of this 
type. 

The idea of the proof is to consider the graph of y = log x and a sequence 
of small sets which lie below y = log x and above the graph of a polygonal 
approximation to y = log x. We show that the areas of these small sets 
form a convergent series and deduce from this that the sequence 


nn+(1/2) 


ün Sr 
O ein!’ 


N= 253500565 
is convergent. Finally we use Theorem F59(b) to find the limit of the 
sequence (an), and this gives Stirling’s Formula. 
Proof of Theorem F62 We need to show that 
n! ~ Varn(n/e)" as n —> oo. 
We begin the proof by considering the function 
f(x) = log z. 


For n = 2,3,..., let P,_1 be the standard partition of the interval [1, n] 
with n — 1 subintervals: 


{[1,2],... [ii +1],..., [n — 1, n]}. 
Now consider the sequence (en) ae where cp is the total area of the 
segments which lie between the graph 

y=logz, TE (1,n), 
and the polygonal graph with vertices 


(1,0), (2,log2), (8,log3), ..., (n,logn), 


4 Stirling’s Formula 


183 


Unit F3 


184 


Integration 


as illustrated in Figure 17. This set consists of n — 1 thin segments. 


®. The function f is concave (that is, its derivative is decreasing), so the 
line segment joining 


(i, logi) to (+ 1, log(t + 1)) 
lies below the graph y = log x. ® 


total area of contribution 
segments = Cy, tO Cn 


(i + 1, log(é + 1)) 


| (i, log i) 


12 Sejid om = 
Figure 17 An approximation for the area under the graph of log x 


The area of the set between y = log x, x € [1,n], and the z-axis is 


/ log z dx = [x log x —x]) 
1 


= nlogn — (n — 1). (23) 
The area between the polygonal graph and the z-axis is 
s(L(f, Pri) +UI Pat) (24) 


®. This is because, on each subinterval [i,i + 1], the function f attains its 
minimum when x = į and its maximum when x = į + 1, so the area of the 
trapezium whose base is the subinterval is equal to the average of these 
two numbers. 8 


Since f is increasing, we have 


L(f, Pr-1) = log1+log24+.---+log(n — 1) 


= log(n — 1)! 
= log(n!/n) 
= logn! — logn (25) 
and 
U(f, Pa-1) = log 2 + log3+---+logn 
= logn!. (26) 


Substituting from equations (25) and (26) into equation (24), we find that 
the area between the polygonal graph and the z-axis is 


5 (log n! — logn + logn!) = logn! — 5 log n. (27) 


4. Stirling’s Formula 


It follows from equations (23) and (27) that 


Cn = nlogn — (n — 1) — logn! + 5 logn 


nrt(1/2) 
= log nn) 
e n: 


PORN . l (1, log 2) (2, log 2) 

The sequence (c,) is positive and increasing. Also, 

Cn <log2, for n € N, 
since the n — 1 segments which contribute to the area c, can be translated 
so that they all lie (without overlapping) in the triangle with vertices 
(1,0), (1,log2) and (2, log 2), as illustrated in Figure 18. This geometric (120) (2,0) 
property holds because the function log is concave. It follows from the i , 
Monotone Convergence Theorem (Theorem D22 in Subsection 5.1 of Figure 18 The segments of 
Unit D2) that the sequence (c,) is convergent. area Cy 


Next, since the exponential function is continuous, the sequence 
ün = Oy = 2, 8,55% 


is convergent also. Thus 


nrt(1/2) 


in = => L asn >œ, (28) 


enin! 
for some non-zero number L. 


To find L, we consider the quotient 


a2 n2n+1 / (2n)2r+(1/2) (2n)! ni/2 


yw a 


ee x 
am €2"—2(nl)2 


e 
enm) (m22 2 
We now let n — œ in this equation. We have 

(2n)! nt? 1 
An > L, am 3 L and (mm a” 
by Theorem F59(b). 


®. Recall that Theorem F59(b) says that 
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Thus we can rewrite statement (28) in the form 
nrn+(1/2) 
MSNu as n > ©. 
Hence, by the Combination Rules for sequences, 
e”n! 
ann 


which can be rearranged to give Stirling’s Formula: 


n! ~ V2rn(n/e)” as n —> ov. = 


=> V27 ano o, 


Summary 


In this unit you have studied the formal definition of what it means for a 
function f to be integrable on a closed interval [a,b]. You did this by 
studying lower and upper Riemann sums for f, which give lower and upper 
estimates for the area between the graph of f and the x-axis fora < x < b, 
by approximating the area with a sum of areas of rectangles. These 
rectangles are based on a collection of subintervals of [a,b] known as a 
partition of [a,b], and you saw that, by decreasing the length of these 
subintervals (that is, the mesh of the partition), we get increasingly 
accurate estimates for the area. If the supremum of the lower Riemann 
sums is equal to the infimum of the upper Riemann sums, then we say that 
the function f is integrable on [a,b] and the common limit is the integral 
of f. You also saw that to prove that f is integrable it is sufficient to 
consider the sequence of standard partitions P, where the interval is 
divided into n subintervals of equal lengths. 


You then met the Fundamental Theorem of Calculus and saw that 
integration can be thought of as the process of finding a primitive, and 
hence as the opposite of differentiation. You studied a number of 
techniques for finding primitives, including integration by substitution, 
integration by parts and the use of a reduction formula. 


You also met some useful inequalities that enable us to estimate various 
integrals when it is not possible to find an exact value. You saw how these 
can be used to prove Wallis’ Formula, an approximation for 7, and studied 
the Integral Test which gives a method based on integration for 
determining the convergence or divergence of certain series. Finally you 
met Stirling’s Formula which gives an approximation for n! and is proved 
using upper and lower Riemann sums. 


Learning outcomes 


Learning outcomes 


After working through this unit, you should be able to: 


e determine the lower Riemann sum L(f,P) and the upper Riemann sum 
U(f, P) for a given function f and partition P 


b 
e understand the definition of the integral f f 
a 
e use upper and lower Riemann sums to determine whether a given 
function is integrable 
e use basic rules for manipulating integrals 
e state various sufficient conditions for a function to be integrable 


e explain what is meant by a primitive of a function and understand the 
Fundamental Theorem of Calculus 


e use the Fundamental Theorem of Calculus and the table of standard 
primitives to evaluate certain integrals 


e use integration by substitution and integration by parts 

e use the reduction of order method to evaluate certain integrals 
e determine lower and upper estimates for given integrals 

e state Wallis’ Formula 


e use the Integral Test to determine the convergence or divergence of 
certain series 


e understand the connection between integration and Stirling’s Formula 
for n! 


e use Stirling’s Formula to determine the behaviour of certain sequences 
involving factorials. 
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Table of standard primitives 


z"+1/(n+1) 
rett Hey + 1) 


Primitive F(x) 


a” /loga 


— COS T 


sin x 


log (sec x) 


e? 


log x 

log |z| 
xlogr—r 
cosh x 


sinh x 


log (cosh x) 


(a? = 22)—1/2, a#0 


(12 — a?)-¥/2, a £0 


(a? + 22) 12, @0 


(a? — 22)1/2, a£0 
(12 — a2)1/2, a£0 
(a? +x)", a A 0 


e% cosbz, a,b #0 


e°” sinbx, a,b £ 0 


{ log(x + (x? — a2)V?) 
cosh !(z/a) 


log(z + (a? + 22)/?) 


tz(a? — 22)V2 + 5a” sin! (x/a) 

5u(u? =a) ža? log(x + (x? — a?) !/?) 
šx(a? +22)V2} 5a" log(a + (a? + x?) 2?) 
a? + b? 


ax 


(acos bx + bsin bx) 


(asin br — bcos bx) 


a? + b2 


Solutions to exercises 


Solution to Exercise F33 
(a) The graph of f is shown below. 


= 0 T 


First, min f = 0, since 

1. f(x) 2 0, for all x € [-1, 1], 

2. f(0)=0. 

Next, inf f = 0, since f has minimum 0 on [—-1, 1]. 

Also, sup f = 1, since 

1. f(x) < 1, for all x € [-1, 1], 

2. if M’ < 1, then M’ is not an upper bound for f 
on [—1,1] because the sequence (1 — 1/n) is 
contained in [—1, 1] and 

f(1—1/n) = (1-1/n)? 41 asn > o, 
so there exists x’ € [—1,1] such that f(a’) > M’. 
(Alternatively, you could consider f(—1 + 1/n).) 


Finally, max f does not exist, since there is no 
point x such that f(x) = 1. 


(b) The graph of f is shown below. 


First, inf f = —1, since 


1. f(x) > —1, for all z € [—1, 1], 


Solutions to exercises 


2. if m’ > —1, then m’ is not a lower bound for f 
on [—1, 1] because the sequence (1/n) is 
contained in [—1, 1] and 

f(1/n) = (1/n)? —1 > -1 asn 00, 
so there exists x’ € [—1,1] such that f(a’) < m’. 

Next, min f does not exist, since there is no point 

x in (—1,1| such that f(x) = —1. 

Also, max f = 1, since 

1. f(x) < 1, for all x € [—1, 1], 

2. f(—1) =1. 

Finally, sup f = 1, since f has maximum 1 

on [—1, 1]. 


Solution to Exercise F34 
The interval [—1, 2] has length 3; hence the 
4 subintervals in this standard partition P of 


[—1, 2] must each have length #. Thus the required 
standard partition of [—1, 2] is 


P ={[-1, -3b l-i b [a4], 12213. 
The mesh of P is the common length of the 
subintervals, namely 3, 


Solution to Exercise F35 


The function and the partition are illustrated 
below. 


eZ 
QOQONONTDS MN) 
SONS 


LLKR 


28 SA NK NK NK NK NK NIKE NIKE 


5 2 


For the three subintervals in P, we have 


mı =0, Mı = f(0) = 1, bx, =}, 
me = f(3)=2, Mo=f(4)=1, 6x =H, 
ms = f (3) =1, M3 = 2, 623 = 5. 


(Note that mg is also equal to f(1).) 
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Then 


3 1 
=la Ta 
_ əl 
~ 50 


Solution to Exercise F36 
Let 


f(z) = zx, TE [0,1], 
and let P, be the standard partition of [0, 1], 


a- (pa pa} 


Since f is increasing we have, for i = 1,2,...,n, 


m=1()= 2, 
“=(= 


and 


OX; = oe 
n 
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and also 


U(f, Pn) = $ Midas 


_ 1 n(n+1) 
NETU 2 


(Here we have used the fact that the sum of the 


arithmetic series 1 +2+---+ n is 


each n € N.) 


Then, as n — oo, we have 


lim L(f, Pa) ae and lim U(f, Pn) 
n— o0 N—+00 


2 


In particular, 


BEREE 


and hence 


f+- e= 


1 
1 
Thus f is integrable on [0,1] and f f= 7 
0 


, for 


z 


Solution to Exercise F37 
(a) The graph of f is shown below. 


YA 

34 e 

0 = 
—2 


Let P,, be the standard partition of [0,1] into n 
equal subintervals: 


(oaae 


For i = 1,2,...,n, we have 
m; =inf{ f(x): (i—1)/n < x <i/n} = —2. 
For i = 1,2,...,n — 1, we have 


M; = sup{ f(x) : (i — 1)/n < x <i/n} = —2. 


Also, Mn = 3 and 


1 
6a, =; fori=1,2,...,n. 
n 


Hence 


Lt) = Som OX; 
i=l 


EC 


i=1 


n—1 


i=l 


5 
=-2+-—>4-2 ano. 
n 


Solutions to exercises 


Since ||P,,|| + 0 as n > oo, and 
lim L(f,P,) = lim U(f, Pa) = —2, 
Noo Noo 


it follows from Theorem F45 that f is integrable 
on [0,1] and 


= 


(b) The graph of f is shown below. 


YA 
H kg 
0 1 2x 


Let P,, be the standard partition of [0,1] into n 
equal subintervals: 


Tesni 


Each subinterval 


li 
E +, for t= 1,2,...,n, 
n n 


contains both rational and irrational points: 


at the rational points, f(x) = zx, 
at the irrational points, f(x) = 0. 


Hence, for i = 1,2,...,n, 
m;=0 and M=— 
n 
Also, 


ih 
6a, ==, fori=1,2,...,n. 
n 
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Hence Solution to Exercise F39 
(f, Pn) = y ee (a) From the Fundamental Theorem of Calculus 
= Nam and the table of standard primitives, we deduce 
n 1 that 
= y (0 x 9 = 0, 4 
i=1 ý f (x? + 9)? dx 
u 0 
4 
U(f, Pa) = >| Mion; = [ea Da a a 
i=1 


= 10 + 2log9 — 3 log 3 


= i 1 9 
=) RA = 10 + 5 log 3. 


i=1 
1 (b) From the Fundamental Theorem of Calculus 
= z2 a and the table of standard primitives, we deduce 
i=1 that 
1 n(n+1) e 
=a 2 f log x dz = [xloga — x]{ 
_i!,1 41 I =(e—e)—(0—1)=1 
= 3 a 6 as n — Oo. : 
Since ||P,|| + 0 as n > oo, but Solution to Exercise F40 
lim L(f,P.) = 07 += lim U(f, Pa), Using the table of standard primitives and the 
ee aces Combination Rules, we obtain the following 
it follows from Theorem F46 that f is not primitives. 


integrable on (0, 1]. (a) F(x) = 4(x log x — x) — 2(5 a) 


Solution to Exercise F38 = 4(xlog x — x) — tan !(x/2) 
e2t 
(a) We have (b) F(x) = 2(} log(sec 3x)) + 52 12 5 (2cos x + sin x) 
1 
F"(x) = ote aie (1 + 3(x? — cr 2) = $ log(sec3x) + že'"(2 cos x + sin x) 
(a? — 4)-1/2 (2? o 4)1/2 4 z) (These results can be checked by differentiation.) 
= 2_ 4yi/2 7 $ 
x + (a? — 4) Solution to Exercise F41 
= (x? — 4)? I 
Fa) (a) Take u = sin 3x; then 
= f(x 
( du 
as reguired. a 3cos3xz, so du = 3cos3z dz. 
(b) We have Hence 
F'(£) = ES cosh x sintim 32) cossadr = + fsmudu 
sinh* x 
= coshx =—icosu 
cosh? z = — + cos(sin 3x). 
= a (b) Taking u = 2 + 3x3, we obtain 
= T 5 
as required. Jre + 32°)" dz = 5.2 + 320). 
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(c) Taking u = 2x7, we obtain 
J esine dx = —} cos(22?). 
(d) Using equation (17), we obtain 


fare + 827) dz = g log(2 + 3722). 


Solution to Exercise F42 
Let u = 1 + e”; then 
du _ 
dz 
Also, 


e”, so du=e’ dz. 


whenx=0, u=2, 
whenx=1, u=1+e. 


Hence 


u e? a du 
——— dr = 2 
» tee w 
u “j 1+e 
"uh 
1 1 e—1 
l+e 2 
Solution to Exercise F43 
(a) Let u = (x— 1)2, so x = u? +1; then 


a =2u, so dx = 2udu. 
du 
Hence 
f dx 
3(x — 1)32 + a(x — 1)1/2 
B f 2u d 
© J 3u + (u? + Iu “ 
= f PE 
4u? +1 
= 2 | _ du 
(2u)? +1 
= tan“! (2u) 
= tan™t(2(x — 1) 2). 


(b) Let u = V1 + e7, so x = log(u? — 1); then 


dr 2u Gu 2u 
a. P ee 


du. 


Solutions to exercises 


Also, 


when zx = 0, u= v2, 
when zx = log3, u=2. 


Hence 
log 3 2 ou 
f evTFEds= | (u? — 1)u 5 du 
0 y2 uł — 1 
2 
-| 2u? du 
V2 
2 
= [28 
= EC im 
= (16 — 4V2)/3. 


Solution to Exercise F44 

(a) Here we use integration by parts, with 
f(e) =logx and g/(x) = 2"; 

then 
f'(x) =1/x and g(x) = 37%, 


Hence 


Ja log z dx = 3 74/3 log r — gf ao dx 


= 34/3 log a — JELO 
= 32% logx— 37%, 


(b) We use integration by parts twice. On each 
occasion we differentiate the power function and 
integrate the trigonometric function. 


We have 
x/2 
f x? cos z dx = Es sin x] an 
0 
T/2 
-f 2x sin x dx 
0 
2 T/2 
a 2 | xsinx dx 
4 0 
and 


x/2 n/2 
f xsin x dx = |x(— cos x)] 
0 


a /2 
-f (— cos x) dx 
0 
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It follows that 


T/2 2 
f t cosrdi = 2 
0 4 


Solution to Exercise F45 


1 
(a) =f e” dx = |e"), = e— 1. 


0 
(b) Using integration by parts, we obtain 


=e—nln-1, forn>1. 


(c) Using the solution to part (b) with 
n = 1,2,3,4 in turn, we obtain 


I =e-Ih=e-(e-1)=1, 
Ig =e —-—2h, =e— 2, 
Iz =e — 3h =e — 3(e — 2) = 6 — 2e, 
Ig = e — 4I; = e — 4(6 — 2e) = 9e — 24. 
Solution to Exercise F46 
(a) Since sin(1/x!?) < 1, we have 
asin(1/x1°) < x, for æ € [1,3]. 


Thus it follows from Inequality Rule (a) that 


3 3 
f xsin(1/x!?) dz < f rdr 
1 1 


(b) If x € [0,5], then 


t= < e?? £ e(1/2)? = en, 


5 Boa a . 
because the function x — e” is increasing on 


10, 51. 


Thus it follows from Ineguality Rule (b) that 


dle 


1m 
< | e” dz < letia 
0 
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Solution to Exercise F47 
(a) Since 

|sin(1/az)| < 1, for x € [1,4], 
and 

2+ cos(1/z) 21, for x € [1,4], 


it follows that 


in(1 
ee <1, forze [1,4]. 
2 + cos(1/z) 
Hence, by the Triangle Ineguality for integrals, 
i 1 
f _sin(1/z) © jsf Jena 
1 2+ cos(1/z) 
(b) Since 


tang > 0, for x € [0,7/4], 
and 

3— sin(x?) >2, for z € [0,7/4], 
it follows that 


tan x 
Saag 125 MAMA 


Hence, by the Triangle Ineguality for integrals and 
Ineguality Rule (a), 


f tang d T tang 
ii H — T 
9 = 3— sin(2?) "|" o = |3— sin(2?) 
T/4 
<f 5 tan xdg 
0 
T/4 


[5 log(sec x)] 5 
= $(log(sec 7/4) — log 1) 
= 1 log(V2) 


a 7 log 2. 


Solutions to exercises 


Solution to Exercise F48 It follows that 
2n)!)? (Qn +1 
(a) We have due gemini? PO EEP 
99 4 22n(n!) 
UK = ms. a a 
2244 64 ((2n)!)?(2n + 1) 
EIUS (2). 
2 24466 256 "\2n+1 
a3 = =. —.-—.—.—- — = — 
3° 13°35 5 T 175° Hence 
b = net 5 > (2n+1 
avr Ng, pe 
294 
bə = (2!)*2 = 4 2, as reguired. 
Al,/2 3 
7 (31)226 16 z Solution to Exercise F49 
3— = zp V9: 
61/3 15 Let u = log x; then 
(b) Using the results of part (a), du 1 
—=-, so du= — 
b? =4 = 3a1, dx x 
9 32 5 Hence 
b3 SG 4342) 
9 2 f dr — fdu 1 1 
p2 256 _ sm x(log x)? u? u log z` 
3 "PSE. , 
75 3 Now let 
as reguired. 1 
(c) We have f(x) = x(log x)? (x € [2, 00)). 
2 (lynn Then f is positive and decreasing on [2, co), and 
no N20." 
VAn n f(x)— 0 as z> o. 
We now try to express 
Also, for n > 2, 
2-2-4-4.-.-- + (2n) (2n) 
Un ee sz HINT n u dx 
[:3.:3-5- = -(n—1)0Qn +1) eS) _ a 
2 2 x(logx)? 
in terms of factorials. 17” 
The numerator is a product of 2n even numbers. = - log -| = 
Taking a factor 2 from each term, we deduce that 1 1 1 
= = £ . 
2:2-4-4..--. - (2n)(2n) log2 logn ~ log2 


Since the sequence of integrals ( | f) is 
2 2 
bounded above, it follows from part (a) of the 


The denominator of a, cannot be treated in quite 
the same way, as all its factors are odd. To relate it 


to factorials, we introduce the missing even factors: Integral Test that the series a 1 > 
e E (Oja (On 1) A. 4 n(log n) 
_1-2-2-3-3-4-4- ++. - (Qn — 1)(2n)(2n)(2n +1) 
= 2-2-4-4.... + (2n)(2n) 
_ ((2n)!)? (2n + 1) 
228 (al)? 
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Unit F3 Integration 


Solution to Exercise F50 


The values are as follows. 


n n! 


30 | 2.652... x 1082 
40 | 8.159... x 10% 
50 | 3.041... x 1094 
60 | 8.320... x 108! 


Solution to Exercise F51 


The values are as follows. 
V2rn (n/e)” 


5 | 118.019... 
10 | 3.598... x 109 
50 | 3.036... x 1094 


Solution to Exercise F52 


In each part we approximate the factorials using 


Stirling’s Formula. 
(a) (300) 1 _ 300! 
150) 2300  150!150! 2300 
6007 (300/e)3°° 


(3007 (150/e)150)? 2300 


v 6007 
3007 


— 30/7 
= 0.046. (to 2 s.f.). 


V6007 (300/e)300 


(100!)3 3300 ~ (2007)3/2 (100/e)300 3300 


= 0.0028 (to 2 s.f.). 
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Solution to Exercise F53 
We have 


Gai = -E 


An / 2n\ _ (4n)!(n!)? 

an} / \n) ~ {2} 
By Stirling’s Formula, and the Product and 
Quotient Rules for ~, we obtain 


(an) / (7) „ Ban (an fey (Vann (r/e) 
2n (Varn (2n/e)2n)° 
B V8rn 44" Irn 
= (4rrn)3/2 26n 
2/8 44" 
~ g 26 
L gn, 


V2 


sO 


n 


Hence A = 1/72. 


Solution to Exercise F54 


Using Stirling’s Formula, we obtain 


n n 


n e 


nm Vrn (nje) yrn 
Thus, by the hint, 


n\ 1/n e” 1/n a 
Gi) EF ( z 7 (Van) /” nin 
We know (from Worked Exercise D26 and 


Exercise D34 in Subsection 3.3 of Unit D2) that, 
for any positive number a, 


a" 51 asn—oo 
and 
n/™ 31 asn— oo. 


Hence 


(vz) +1 as i oo 


and 


ni/n +1 asn 3 oo. 


Solutions to exercises 


It follows that 


n” 1/n 
( ) ~e an oO; 


n! 


that is, 


197 


